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The detailed derivation of a theory for transport in quasi-two-dimensional metals, with small- 
angle elastic scattering and angle-independent inelastic scattering is presented. The transport equa- 
tion is solved for a model Fermi surface representing a typical cuprate superconductor. Using the 
small-angle elastic and the inelastic scattering rates deduced from angle-resolved photoemission 
experiments, good quantitiative agreement with the observed anomalous temperature dependence 
of the Hall angle in optimally doped cuprates is obtained, while the resistivity remains linear in 
temperature. The theory is also extended to the frequency-dependent complex Hall angle. 

PACS numbers: 74.20.-z, 74.20.Mn, 74.25.Fy 



I. INTRODUCTION 

Soon after the discovery of high- temperature supercon- 
ductivity in the cuprate compounds, it was found that all 
the transport properties for compositions near those for 
the maximum T c have anomalous temperature and/or 
frequency dependence, in contrast to what is expected 
for Landau Fermi liquids. Most of the anomalies, in re- 
sistivity and optical conductivity, in the frequency de- 
pendence of Raman intensity, in the tunneling spectra, 
and in the temperature dependence of the copper nu- 
clear relaxation could be understood by the "Marginal 
Fermi Liquid" (MFL) phenomenology^ Within MFL, 
scale-invariant fluctuations in both in the charge and 
magnetic sectors are assumed to have the form, 



Imx(q,w,T) 



-N (w/T), 
-7Vo(sgnw), 



cj<T 

T<C^<Cw c , (1) 



characteristic of a quantum critical point, uj is the fre- 
quency and q the momentum of the fluctuations, Nq 
is the density of energy states per unit volume and lu c 
is a high-frequency cutoff. At long wavelengths this 
form is observed directly in Raman scattering^. A prin- 
cipal prediction^ from this hypothesis is that at low 
energies (u> <C T), the inelastic part of the single- 
particle relaxation rate has the MFL form AT with co- 
efficient A having negligible momentum dependence ei- 
ther along or perpendicular to the Fermi surface. This 
form has been confirmed in angle-resolved photoemission 
(ARPES) experiments^ and leads directly to the observed 
linear temperature dependence of the resistivity. 

However, the temperature dependence of anomalies in 
the normal state Hall effect and magnetoresistance could 
not be understood by the MFL hypothesis alone. For ex- 
ample, for a situation near optimal doping, where the re- 
sistivity is linear in temperature, the expectation is that 
the cotangent of the Hall angle (a xx /a xy ) should also be 
linear. Experiment shows it to be more nearly quadratic. 
This left open the possibility that essential new physics 



near optimum doping may not be captured by the MFL 
scenario. In this paper we show that the temperature 
dependence of the Hall angle may be understood quanti- 
tatively by a proper application of transport theory using 
the measured single-particle relaxation rates. 

The single-particle self energy E(k, co,T) is mea- 
surable in angle-resolved photoemission (ARPES) 
experiments^ As stated above, within MFL, the predic- 
tion is that the inelastic part of E is k-independent and of 
the MFL form, 1 proportional to ui for u 3> T and to T for 
T ^S> uj. ARPES experiments^ do find the inelastic part 
E of this form but find in addition an elastic part which 
varies in magnitude around the Fermi surface. Thus, as 
explained in detail earlier^ on the Fermi surface (u> = 0), 
the experimentally-measured self energy consists of the 
MFL part AT and an anisotropic T-independent elastic 
part: 



ImE(fc,T) = AT + 7 (fc), 



(2) 



According to experiment, j(k) increases by about a 
factor 4 to 5 going from the (tt, tt) to the (tt, 0) direc- 
tion along the Fermi surfaced A crucial point is that 
even at its minimum value, it is more than an order of 
magnitude larger than the transport scattering rate due 
to impurities obtained by extrapolating the normal state 
resistivity to T — 0. As argued in Ref. 4, such behavior 
arises if j(k) comes from small-angle impurity scattering. 
We discuss this point further in Sec. VII. 

In a previous communicationj 5 we described, for high- 
Tc superconductors, how Eq. (I) can be used in a 
Boltzmann equation analysis to account for observed 
anomalies 6 7 8 910-11 in the Hall effect. We performed a 
calculation (page 4655 of Ref. 1) using a simple Fermi 
surface in order to give an example of how a new 
contribution's could dominate the conventional result 
and thus account for experimental observations. However 
it was pointed out to us by V. Yakovenkoi^ that we erred 
in the form we chose to parameterize j(k)£m The pur- 
poses of this paper are to give a more complete solution 
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of the Boltzmann equation, correct the above error, and 
to present a model calculation which illustrates how, with 
a proper parameterization and a reasonable choice of pa- 
rameters, the presence of the anisotropic impurity scat- 
tering together with other anisotropics can account for 
the observed behavior of the Hall angle. 6- '• 8 - 9 - 10 - 11 po- 
ther, we generalize our results to the complex Hall angle 
at finite frequencies which has been recently measured, 15 

A point about the experimental results needs stressing. 
In many presentations of the data, for example Refs. 9,10, 
single temperature non-integer power-law fits to the Hall 
number or the Hall angle have been made, with vary- 
ing success, while integer power law fits£ do not usually 
faithfully represent the data. There is no physical reason 
to expect integer power law behavior and our analysis 
shows quite generally that the solution of the Boltzmann 
equation can give a sum of contributions with different 
temperature dependences. In particular, it is usually 
found2*iii that a temperature power law with a power 
less than 2 can fit the data for the cotangent of the Hall 
angle. As we show here, this is indistinguishable from 
the ratio of terms each of which is a sum of contributions 
with different temperature dependences. 

Our results have two significant conclusions. First, 
anisotropies, especially in the scattering rate, introduce 
corrections to the calculation of magnetotransport prop- 
erties so that conventional results such as Rh = (nee) -1 , 
or tan 9h = u c T tr are not in general valid>i& Second, the 
fact that the Hall effect in the cuprates can be understood 
within transport theory, with only the measured single- 
particle relaxation rate, the measured resistivity and 
Fermi surface quantities as inputs implies that no new 
physics is involved other than what leads to the behavior 
of those quantities. We conclude that the MFL, which 
gives the experimentally observed temperature depen- 
dence of the relaxation rate, is sufficient to account for all 
the the anomalous magnetotransport near optimal dop- 
ing. This point of view is reinforced by recent frequency- 
dependent complex Hall effect measurements^ whose re- 
sults are consistent with the analysis presented here, as 
discussed in Sec. VI. 



II. TRANSPORT EQUATION 

At long wavelengths and low frequencies the deriva- 
tion of Boltzmann equation relies only on conservation 
laws and is valid whether the system is a Fermi liquid or 
not. Therefore we start with the linearized Boltzmann 
equation for the deviation <7(k, t) of the momentum dis- 
tribution function from its equilibrium value /(k) in the 
presence of uniform static electric and magnetic fields, 



— + eE v k — 

dt <9e k 



h c 



(v k xB)-V k .g k = C k (3) 



Here C(k) is the collision operator. The solution to this 
equation for the stationary case (uj — 0) is given bj*i£ 



k,k' 



E • v k ' ~ 



df_ 
<9ek' 



where 



j4 k ,k' = h 



7^) + ^ VkXB ' Vk 



<5 k , k ' — C] 



k,k' 



(4) 



(5) 



Here C k . k ' is the "scattering-in" term in the collision 
operator for the Boltzmann equation and 7i/r(k) = 
X] k / Ck,k', is the "scattering-out" term equal to the 
single-particle relaxation rate. It is evident that the dis- 
tribution <?(k) is not only determined by the energy of 
the state k but also by the anisotropy of the scatter- 
ing. The distribution is depleted in directions of large 
net scattering and augmented in directions of small net 
scattering. 

We calculate the conductivities using the single- 
particle scattering rate of Eq. (2). The conductivity ten- 



sor is 



t'" 7 



n 



E 

k,k' 



Vk' 



(6) 



where fl is the sample (taken to be a plane) area. We 
expand A -1 from Eq. (4) in powers of B>i£ 

A~ 1 =T - (e/c)T(v x B • V) T 

+ (e/c) 2 T(v x B • V) T (v x B • V) T, (7) 

where (from Eq. (2.12) of Ref. 17) 

7k, k' = T-[r k Jfc.k' + ^2 r k Ck.k" 7k" ,k'] (8) 
11 k" 

The first term in A^ 1 gives the longitudinal conductivity, 
the second the Hall conductivity a xy and the magneto- 
conductivity may be calculated from the third term. 

As discussed above, the r k is the MFL scattering plus 
the angle dependent impurity piece: 



l/r k = 1/tm + l/n (k). 



(9) 



Our approach to solving Eq. (7) is to take advantage 
of the properties of 1/tm and 1/tj. The kernel C in 
Eq. (8) comprises the vertex corrections for the various 
scattering mechanisms. Within MFL, the 1/tm is an s- 
wave scattering process so the MFL interaction does not 
appear in C. That leaves 



C k ,i 



2^(e k -e k 0MM')| 2 



(10) 



where, e.g., 9 is the angle of k along the Fermi surface. 
We have argued^ that the impurity scattering matrix 
element Ui(9, 9') involves scattering through small angles 
only. This property allows an expansion of the RHS of 
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Eq. (8) in powers of a "small angle scattering parameter" , 
9 C which can itself be a function of 9. We illustrate such 
an expansion for the calculation of 1/rf. 



-(0) = ~£k(M')| 2 <Ke'-e) 



M(9,9)N(9) + ^ 



Here we have used 



n dk t /de 

h v{6) ' 



(11) 



(12) 



N(9) is the density of states per unit energy at the Fermi 
surface at angle 9 and v{9) is the Fermi velocity. dk t 
is an infinitesimal taken tangent to the Fermi surface at 
the angle 9. We have abbreviated \u,{9,9')\ 2 by U(9,9'). 
We shall describe below how experiments show that 9 C is 
indeed sufficiently small to justify such an expansion. 



III. LONGITUDINAL CONDUCTIVITY 

We begin by determining the longitudinal conductivity 
within the model. We first give the solution for a general 
Fermi surface and then specialize to a particular example 
to illustrate the results. 

From Eqs. (5,6), we have 



kk' 



n 



Define the vector L as 

L k =^T(k,k')(-d//de')v k ,, 



so that 



xx c \ x T x 



(13) 



(14) 



(15) 



Using Eq. (8) in Eq. (14), we find the integral equation 
for L: 



T - Tk 



(-T^K+E^k"!*" 



de 



(16) 



As discussed earlier, C contains only the small angle scat- 
tering so that 



T - Tk 



(-^W + / d9'N(9')U(9,9')L h 



(17) 



which, since L is restricted to the Fermi surface, can be 
rewritten as 



+ 



d9'N(9')U{9,9')[L(9') - L(0)], (18) 



where 1/r, = 1/r — 1/tm is given in Eq. (11). Given 
the small-angle scattering restriction on U(9,9'), we ex- 
pand the difference L(0') — L(9) and find the differential 
equation (primes indicate derivatives with respect to 9) 



M(9) = v(0) + u(9)M"(0) + 2u'(0)M'(0), 



where 



M(9) = 



r M (-df/de 



MO), 



(19) 



(20) 



and 



u(9) 
u'{9) 



(9 3 c /24)(r M /h)N(9)U(9,9), 
(9l/24)(r M /h)[^N(9')U(9,9%^ 9 .(21) 



This differential equation can be solved for the compo- 
nents of M, hence L, once the 0-dependences of v, u and 
u' are known. This equation, which is basic to the evalu- 
ation of both the longitudinal and the Hall conductivities 
is completely equivalent to Eq. (10) of Ref. 5. Although 
u contains a small parameter an iterative solution of Eq. 
(19), as in our previous workj£ is not valid. We solve Eq. 
(19) exactly in this paper. 

As we discussed previously^*i^ the details of the 
anisotropies of the quantities entering the transport co- 
efficients determine the magnitudes of the various contri- 
butions, in particular to a xy in the presence of a magnetic 
field. The shape of the Fermi surface is especially impor- 
tant since it determines the size of the contribution to 
a xy from the isotropic marginal Fermi liquid scattering 
rate 1/tm- In what follows, we take the anisotropy of the 
impurity scattering from ARPES data and we assume a 
simple form for the Fermi surface velocity. 

A schematic of the Fermi surface is shown in Fig. 1, 
where the angular variable 9 is shown. A general Fermi 
surface respecting the square symmetry of the Cu02 lay- 
ers must have Fermi velocities in the first quadrant of the 
form 

v x = sin[(2n+ 1)9] 

n 

v y = ^ n (-l) n cos[(2n+l)0]. (22) 

n 

For our calculations we take a form for v(0) which is 
the simplest extension beyond the circular Fermi surface 
(similar to that shown in Fig. 1) consistent with Eq. (22). 
In the first quadrant, 



uc^sin^ + psin 39) 
vq(cos 9 — p cos 39) 



(23) 



For this choice for v(9), the density of states of Eq. 
(12) in all quadrants is 



2nh v (l - /5C0S46)) 5 / 4 ' 



(24) 
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FIG. 1: Schematic Fermi surface. 

where fco is the value of &f(6>) at 9 = 0. 

Thus, according to Eq. (15), a xx is given by: 



e 2 r M 



(1-p) 1 / 4 fc f^ 2 d() _ v x (9)M x (9) 



ko r 

vo Jo 



(1 -pcos46») 5 /4- 



(25) 



IV. HALL CONDUCTIVITY 



For a magnetic field perpendicular to the xy plane, 
v x B • V = [fl/2irhN(9)}B ■ d/d9 and from Eqs. (6,7), 



2^c~ ^ 



k.k'.kl 



A^l) 



We rewrite this using Eq. (14): 



e 3 B 
2^c~ 



k,kl 



T(k,ki) d 

AT(6»i) d6i 



We define 



so that 



p 3 

_ ^ \ x K y 



(26) 
(27) 

(28) 

(29) 



Using Eq. (16) in Eq. (31), we get an integral equation 
for K: 



= |(4w + E^K k „| 



L k ' 



h \n(0) 



d9'N(8')U(6,6')K(6') } .(30) 



We proceed exactly as in the analysis for L, Eqs. (16- 
21). For the required K y , we find the differential equation 

Z{9) = M' y {9){\ - pcos40) 5 / 4 + u(9)Z"(9) + 2u'{9)Z'{9) 

(31) 

where Z(9) determines K y (9) as 

and My is obtained as the solution of Eq. (19). As in 
that case, an iterative solution of Eq. (31) is not valid. 
Combining all our results into Eq. (29), we find 



tt/2 



d9 



e 3 B f T M 

' tt 2 c\ h J J (1 -pcos46») 5 / 4 



v x {6)Z{6) 
(33) 



When there is no small-angle scattering, M y — v y and 
Z(8) = -v (sm8 - 3psin30)(l - pcos46») 5 / 4 . In that 
case, 



,7T/2 

a xy oc / d9(sin 2 - 3p 2 sin 2 9) 
Jo 



(34) 



which vanishes when p = 1/ >/3. This is the familiar re- 
sult of the vanishing of a xy when the Fermi surface con- 
sists of equal portions of positive and negative curvature. 
One sees that the anisotropic small angle contributions, 
which have a leading contribution proportional to r^, 
i.e. to 1/T 3 , do not vanish there. 



V. CHOICE OF PARAMETERS AND 
COMPARISON WITH EXPERIMENTS 

We now describe the evaluation of the conductivities 
using experimental data from ARPES and from the lon- 
gitudinal transport in zero magnetic field. The strategy 
is as follows. The anisotropics in the problem will be de- 
termined from the anisotropy of the impurity scattering 
l/n(9) = 2j(k)/h as determined from ARPES. For sim- 
plicity, we shall let the density of states N(9) be responsi- 
ble for the anisotropy of 1/rj, and hence for the quantity 
u(9), see Eq. (21). The measured anisotropy of 1/rj from 
the 7r, 7r direction to the 7r, direction in the Brillouin 
zone then determines the velocity anisotropy parameter 
p from Eqs. (11,24). These parameterizations give the 
correct behavior at the edges of the Brillouin zone. The 
small angle parameter 9 C can then be found from the 
"residual resistance" ratio (RR) of the resistivity at T c 
to the extrapolated value at T = 0. That is sufficient to 
determine the T-dependence of cot6># = a xx /a xy . The 
magnitude of cot 9h is then fixed by the effective mass 
hko/vo, which we leave as an adjustable parameter. It 
will be seen that the effective mass has a quite reasonable 
value. 
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Resistivity 



50 100 150 200 250 300 350 



FIG. 2: Linear temperature fit to calculated resistivity (dots) 
from Eq. (33) for 2 C /2A = 0.0014. The resistance ratio, 
p(100)/p(0), is 8. 



The conductivities are given by Eqs. (25,33). The 
quantity u which appears in the differential equations 
for M and Z is defined as u(9;T) = {6 2 C /2A)T M /Ti{0) = 
(tm /h){Ql/2A)N{0)U{6,0). The measured anisotropyi 
of 7(fc) = h/2ri is about 4.5. We assign this to N(9), 
thus determining p from Eq . (24) as p = 0.54. There- 
fore, we have 



eta 

O 

o 
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1/ Ti (9) = (l/r o )(l-pcos40)- 5 / 4 , 
u(0;T) = w (T)(l-pcos46»)~ 5 / 4 , 
p = 0.54. 



FIG. 3: Data of Chien, et al, Phys. Rev. Lett. 67, 2088 (1991). 



(35) ous result for a xx (T) to give cot Oh (T) = o xx (T) j o xy (T): 



The quantity u (T) = r M (T)(0 c 2 /24)/r o . r M (T) and r 
may be determined from ARPES as described in Ref. 3. 
From ARPES, h/n w 0.24 eV at 6 = tt/8. So ?j/t = 0.24 
eV. From the same source, %/tm = 0.015(T/100) eV, 
where T is in Kelvin. This determines uq as uq{T) = 
16(lOO/T)(0 2 /24). The remaining parameter uq, that is 
9 Cl is determined from the resistance ratio RR as follows. 

The longitudinal resistivity in zero magnetic field p xx 
is given by the inverse of o xx from Eq. (25). Unfortu- 
nately we do not have an analytic solution of the dif- 
ferential equation Eq. (19) for M x , although the coeffi- 
cients are now known from Eq. (35). Numerical integra- 
tions of Eqs. (19,25) are performed for different values of 
Uq(T) = (1600/T) (02/24). For a given value of 6 C , p xx is 
almost precisely linear in T above T = 100 K. See Fig. 2. 
The T = intercept of p xx , po, is obtained by extrapolat- 
ing the found high temperature linear dependence. The 
RR is defined as p xx (T — 100)/poj and is typically about 
8, say. This is obtained for u (T) « .022(100/T) so that 
the characteristic small angle parameter is indeed small, 
6*2/24 ps 0.0014. 

All parameters are now fixed with the exception of 
to* = hko/vo. The numerical integrations of Eqs. 
(19,31,33) for M y , Z and a xy respectively can be carried 
out for different values of T and combined with the previ- 



cot 9 H {T) 



0.82 to* J d9n(9)v x (9)M x (9) 



io c TM m J d9n(9)v x (9)Z(9) 



(36) 



where we have used p = 0.54, oj c is the cyclotron 
frequency with the bare mass and n{&) = (1 — 
0.54 cos46»)- 5 / 4 . 

We now compare our result with the experiment of 
Chien, et al£ They measured the Hall angle as a function 
of temperature for various concentrations of Zn impuri- 
ties in near-optimally doped YBCO. Their data at B = 8 
T (u c = 1.28 x 10 12 s" 1 ) is reproduced in Fig. 3. 

We examine the data for zero Zn concentration. In 
Fig. 4, we have replotted the x = points from Fig. 
3 (large dots) and the theoretical result (small dots) of 
Eq. (36) with m*/m = 1.5, so adjusted to give the best 
agreement. It is important that to* turns out to be rea- 
sonable. We emphasize that the parameters we picked 
were determined from ARPES experiments on BiSCO, 
while the Hall data is on YBCO. Thus, our intent here 
is only to show that the scattering rates characteristic of 
high-T c superconductors in the normal state near opti- 
mal doping can by themselves account for the temper- 
ature dependence of the Hall angle, not to demonstrate 
quantitative agreement with a particular experiment. It 
is seen that the experimental temperature dependence is 
fairly well reproduced. 
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FIG. 4: Theory fit to x = data of Chien, et al£ 
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FIG. 5: Comparison of theory (small dots) including in-plane 
scatterers with data (large dots) from Ref. 9, for optimally 
doped BSLCO. 



Notice that the experimental curves have a negative 
curvature on a T 2 plot. While this might suggest a power 
law less than 2 and several investigators have tried such 
a fitjSiiS our picture indicates that this behavior is a con- 
sequence of the fact that the conductivities are each the 
sum of two terms with approximate l/T™ and \/T n+1 
(a xx : n — 1, a xy : n = 2) dependences. 

The data of Ref. 6 and our theory both suggest that 
in the absence of in-plane impurities such as Zn, the ex- 
trapolated cot Oh at T — > is zero. In-plane impurities 
will in general not give scattering restricted to small an- 
gles. Indeed, their influence is evident in the measured 
resistivity. We could include their effects by assuming 
them to be isotropic scatterers which add a constant 1/t z 
to the isotropic linear in temperature MFL rate 1/tm- 
This has the effect of an almost parallel upward shift of 
cot 9h so that it has a zero temperature intercept pro- 
portional to 1/t z , i.e. to the in-plane impurity concen- 
tration. We illustrate this qualitatively using some data 
from Ref. 9. In Fig. 5, we show the data (large dots) 
for Bi2Sr 2 _ s La z: Cu06 with x — 0.44 (optimal doping). 
It is seen that the data does not extrapolate with a zero 
intercept. We interpret this as signalling the presence of 
in-plane scatterers. The theoretical curve (small dots) 
uses similar parameters as in Fig. 4 but with an in-plane 
impurity contribution added to 1/tm- 



VI. COMPLEX HALL CONDUCTIVITY 

In this section, we comment on the recent ac Hall ef- 
fect results of Grayson, et al^ For low frequencies and 
T > T c , we can replace Itm ~ 1/ T tr by 1/ T ti' — This 
follows immediately from Eqs. (3-5): the distribution 
function g(k., ui) and therefore the transport properties 
are obtained from the w = results by this replacement. 

From Fig. 2 it is seen that the resistivity is almost 
precisely linear in T in the normal state: p xx = a + bT, 
where a is the (small) residual resistivity due to the small- 



angle impurity scattering. Thus, a/6 is about 10, say. 
We have also found that a xy (T) is, in the temperature 
range of interest, very closely of the form c/T 2 + d/T 3 , 
with c/d « 0.01 K _1 . Therefore, neglecting quantities of 
order 0.1, our prediction for the Hall angle is 



tan(9 ff = a xv p xx 



A 



1 



1 



(1/T tr - iuj) 2 d(l/Tt r — Uj) 



(37) 



where A and c/d are constant in temperature. The latter 
ratio is now in units of seconds; for our example param- 
eters it is about 2 x 10~ 14 s. We emphasize that a single 
relaxation rate l/r tr enters all our expressions. The con- 
stant A in Eq. (37) corresponds to the quantity ujh^p 
of Ref. 15. In that reference, Fig. 4 shows that A is in- 
deed temperature and frequency independent. The c/d 
term is the conventional term. In Ref. 15, it was shown 
that by itself it cannot account for the experimental data. 
Just as the dc data shows deviations from T 2 behavior 
for cot#jf, we expect that the ac data should be fit by a 
combination of the two terms in Eq. (37) . At higher fre- 
quencies or temperatures, the conventional term should 
finally dominate. The complex Hall angle measurements 
convincingly demonstrate that just one ineleastic trans- 
port rate determines all the frequency and temperature 
dependence of the transport properties of the high-T c su- 
perconductors in the marginal Fermi liquid region, i.e. in 
their normal state near optimum doping. 

When calculating the ac conductivity, the limits of ap- 
plicability of the Boltzmann equation should be kept in 
mind. The results are only valid for lot < 1, In this 
regime the conservation laws (in terms of bare particles) 
completely determine transport and the Boltzmann equa- 
tion deals with them properly. In a microscopic theory, 
the effects of small-angle scattering calculated here ap- 
pear as corrections to the coupling of the external mag- 
netic field to the carriers. At high enough frequency these 
vertex corrections must vanish. The crossover frequency 
and the behavior of transport properties in the interme- 
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diate regime can only be determined by a microscopic 
calculation. Experimental results at high frequencies do 
depart from the predictions of Eq. (37). 

VII. CONCLUDING REMARKS 

Here we discuss sundry issues related to the theory 
and calculations presented in this paper. We emphasize 
that this work pertains only to the Hall effect in opti- 
mally or overdoped samples. For the latter, the linear- 
in- temperature scattering rate 1/tm should be replaced 
in the theory by a rate with the T-dependence of the 
observed resistivity. In undcrdoped cuprates, additional 
physical considerations appear to determine the trans- 
port and equilibrium properties. 

A. Magnetoresistance 

Earlier, 5 we gave a plausibility argument that the 
observed magnetoresistance would follow quantitatively 
from the same considerations which lead to the observed 
Hall effect. We leave it as a future exercise to calculate 
this quantity directly from the solution of the Boltzmann 
equation as an extension of the present theory. 

B. Small-angle scattering, Fermi surface geometry, 
bilayer splitting . . . 

We have modeled the Fermi velocity through Eq. (23), 
which with p w 0.54 is in itself enough to give the ob- 
served angular variation of the elastic part of the single- 
particle relaxation rate. As we have seen, this also gives 
agreement with the observed temperature dependence of 
the Hall angle. This value of p is 10% away from the 
value at which the customary Hall coefficient would be 
zero. This seems to us reasonable, since in a model with 
nearest neighbor Cu-0 hopping alone, the Hall effect is 
zero at half-filling and with 0-0 hopping included, the 
density for zero Hall effect shifts substantially towards 
hole doping. The optimal doping composition lies in the 
range 15-20% hole doping. In this connection, it is impor- 
tant to remember that we are using parameters obtained 
from experiments on BSCCO to fit Hall data on YBCO. 
Therefore the precise numbers are not too meaningful; 
we only argue that we have shown that the small-angle 
scattering scenario gives the qualitative features of the 
experiments. 

In our numerical analysis, we have chosen the simplest 
parameterization of the Fermi velocity to show the plau- 
sibility of the small-angle effect in relation to the experi- 
ments. Our parameterization should not be expected to 
represent the Fermi surface adequately. For the actual 
Fermi surface several terms in Eq. (22) for the Fermi ve- 
locity would have to be included with coefficients smaller 
than p. We note however that the successive terms give 



compensatingly larger coefficients since the higher har- 
monics produce larger derivative terms in Eq. (36) for 
a xy . For a truer Fermi surface the advantage of analyt- 
ical calculations is lost. These remarks are of academic 
interest at the present time since the actual Fermi surface 
is not accurately known. 

More relevant is the fact that the Fermi surface of Bi- 
2212, the ARPES data for which we have used here, has 
two sheets coming from the bilayer splitting. The bilayer 
splitting varies as a function of angle, being largest in the 
7r, direction. This arises because of the geometry of the 
interlayer orbitals4& This angular variation is essentially 
the same as that of the observed anisotropic elastic scat- 
tering. The fact that the bilayer splitting is not resolved 
in the experimental data we have used indicates the pres- 
ence of an elastic scattering mechanism which couples the 
layers. We surmise that for interlayer or intralayer scat- 
tering due to impurities between the planes, the same 
orbitals are involved and therefore the same angular de- 
pendence arises in scattering. In a single-layer model 
as treated here, this effect is parameterized through the 
choice of the angle-dependence of the Fermi velocity as 
in Eq. (23). 



C. Comments on related work 

As already remarked, although the analytical expres- 
sions derived in our earlier work£ are correct, we erred, as 
pointed out by V. Yakovenko, in our choice of the param- 
eterization in which they were evaluated. This error has 
been remedied here. Our more complete analysis with a 
simpler parameterization in fact gives results consistent 
with experiment. We have shown that not only do we get 
the right temperature dependence but that with a reason- 
able effective mass, the absolute value of the cotangent 
of the Hall angle is obtained within 10% of the data from 
the same parameterization. 

In Ref. 19, R. Hlubina argues that the new term discov- 
ered by us, if adequate to explain the Hall angle, would 
make an unacceptable contribution to the longitudinal 
resistivity. Since our parameters are actually determined 
by the measured resistivity, this criticism is invalid. It 
is difficult to make a direct comparison with the calcula- 
tions of Ref. 19 since a quite different parameterization 
of the Fermi surface is used there. However, we can note 
that the parameterization family used in Hlubina's eval- 
uation does not include ours and gives a much larger 
conventional contribution to a xy than the one we have 
used. 

A recent paper of Carter and Schofield 20 addresses the 
main point of our original paper£ Their work consists of 
two parts. One is analytical, the other is a numerical so- 
lution of equations we derived in Ref. 5. We cannot com- 
ment on the adequacy of the numerical work, but Carter 
and Schofield arc indeed correct that the customary term 
(with a different temperature dependence than the new 
term) is zero only for a particular choice of Fermi surface. 
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This is of course well-known: the traditional Hall angle is 
zero only for "particle-hole symmetry" defined in terms 
of the curvature of the Fermi surface. However, as we 
have shown here, it is not necessary that the customary 
term be zero (and all the temperature dependence come 
from our new term) to get good agreement with experi- 
ment. It is sufficient that the customary term is about a 
factor of four smaller than the new term at T w lOOiV 
In fact, the experimental data for cot Or, when plotted 
against T 2 , usually shows a slight downward curvature. 
This point has been discussed in Sec. I and is adequately 
demonstrated in the experimental data shown in Fig3. 
3-5. 

The analytical calculation of Ref. 20, expressed in their 
Eqs. (7-10), is not sufficiently general to contain the 
anisotropies required to show the effect we have derived. 
Their calculation for a xy is equivalent to retaining only 
the impurity scattering contribution to the conventional 



term, not the new contribution which we have derived. 
See footnote 22 of Ref. 20 for a comment on this. Ac- 
tually, the analytical calculation of Ref. 20 is based on a 
circular Fermi surface and we agree that such a parame- 
terization never leads to a large enough effect. 
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